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SECTION  1 


INTRODUCTION 


This  Research  and  Development  Technical  Report  covers  work  per- 
formed by  The  Charles  Stark  Draper  Laboratory,  Inc.,  (CSDL)  for  the 
U.S.  Air  Force  Systems  Command,  Aeronautical  Systems  Division,  Wright- 
Patterson  Air  Force  Base  under  Contract  F33615-77-R-1103  between 
1 February  and  31  December  1977. 

The  effort  performed  under  this  contract  is  a natural  follow-on 
to  a CSDL  effort  conducted  during  calendar  year  1976.  Reference  1 
documents  the  work  performed  under  this  initial  effort.  Its  objective 
was  to  develop  the  tools  and  techniques  required  to  assess  the  techni- 
cal feasibility  of  using  a minimum  number  of  inertial  sensors  to  pro- 
vide the  inertial-reference  information  consistent  with  the  weapon- 
delivery,  flight-control,  navigation,  fire-control,  and  flight  safety 
requirements  of  an  aircraft.  The  term  Multifunction  Inertial  Reference 
Assembly  (MIRA)  describes  the  generic  class  of  systems  designed  to 
meet  this  objective. 

The  initial  MIRA  effort  was  primarily  oriented  toward  developing 
and  demonstrating  simulation  capabilities  to  evaluate  MIRA-type  systems 
in  realistic  environments,  and  toward  considering  the  life-cycle-cost 
implications  of  the  MIRA  approach.  The  main  goal  was  to  develop  a ver- 
satile evaluation  capability  that  would  have  long-term  utility  in  the 
MIRA  effort. 

The  present  CSDL  effort  documented  in  this  report  is  denoted 
Multifunction  Inertial  Reference  Assembly  Technology  (MIRAT)  to  differ- 
entiate it  from  the  first  effort.  It  has  three  main  objectives.  The 
first  is  a continuation  of  the  simulation  development  begun  under  the 
first  CSDL  MIRA  effort.  Section  2 of  this  report  documents  the 
results.  The  tasks  performed  during  this  phase  of  the  contract  were 
the  development  of  a real-time  digital  simulation  with  the  same 
capabilities  as  the  nonreal-time  MIRA  simulation  developed  by  CSDL 
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under  the  first  MIRA-related  effort,  and  the  addition  of  an  air-data 
model  and  a transport-vehicle  model  into  the  simulation.  Also  incor- 
porated was  a navigation  and  attitude  error-propagation  scheme  using 
incremental  velocity-  and  angular-sensor  errors  as  inputs  together 
with  flight-profile  data.  This  algorithm  makes  it  possible  to  isolate 
the  effects  of  sensor  errors  on  navigation  performance.  Two-degree- 
of-freedom-gyro  models  were  also  added  to  the  MIRA  simulation  to  comple- 
ment the  existing  laser-gyro  models. 

The  second  major  CSDL  task  performed  under  this  contract  was 
directed  toward  the  utilization  of  the  simulation  tools  in  the  overall 
MIRA  test  and  evaluation  program.  Parametric  simulation  studies  were 
performed  to  evaluate  the  effects  of  sensor  location,  sensor  type, 
sensor  errors,  etc.  upon  the  navigation  and  flight-control  capability 
of  the  aircraft.  In  addition,  further  simulation  studies  were  performed 
to  evaluate  specific  MIRA  configurations  and  provide  an  independent 
assessment  of  candidate  configurations.  The  configurations  to  be  evalu- 
ated were  to  be  given  to  CSDL  based  upon  the  results  of  the  McDonnell 
Douglas  Aircraft  Company  MIRA  effort.  However,  this  information  was  not 
supplied  to  CSDL  in  time  to  be  taken  into  consideration  on  the  present 
program.  Therefore,  CSDL  selected  the  configurations  to  be  evaluated. 
Section  3 describes  the  simulation  runs  performed,  and  presents  the 
results  and  conclusions  obtained  from  them.  The  evaluation  of  the  MIRA 
configurations  selected  by  the  McDonnell  Douglas  Aircraft  Company  will 
be  completed  under  an  extension  to  the  present  contract.  The  results 
of  this  evaluation  will  be  published  in  an  appendix  to  this  report. 

The  last  major  task  required  for  CSDL  under  this  contract  was 
the  development  of  a test  and  evaluation  plan  for  MIRA-type  systems. 

CSDL ' s proposed  plan  is  presented  in  Section  4. 

Section  5 summarizes  the  total  CSDL  MIRA-related  effort,  and  pro- 
vides information  on  future  MIRA  use  of  the  tools  developed  and  results 
obtained. 
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SECTION  2 


MIRA  SIMULATION  DEVELOPMENT 


One  of  CSDL's  major  goals  on  the  MIRA  program  was  the  development 
of  a simulation  capability  to  evaluate  the  MIRA  concept.  It  was  de- 
cided that  the  simulation  should  be  completely  digital.  Furthermore, 
both  nonreal-time  and  real-time  versions  were  developed.  The  nonreal- 
time version,  coded  in  FORTRAN,  is  simpler  to  implement  and  allows  for 
easier  checkout,  while  the  real-time  version,  coded  in  assembly  language, 
permits  human  and/or  equipment  interaction  with  the  simulated  system. 

The  initial  CSDL  MIRA  effort  concentrated  on  the  development  of  the 
nonreal-time  simulation  capability,  while  the  development  of  a parallel 
real-time  capability  was  emphasized  during  this  contractual  effort. 

Another  primary  feature  of  the  MIRA  simulation  is  flexibility.  Provi- 
sions have  been  included  in  the  simulation  to  allow  parameters,  sensor 
configurations,  sensor  locations,  program  subroutines,  models,  etc.  to 
be  changed  easily. 

The  basic  computation  unit  was  a Xerox  Data  Systems  (XDS)  9300 
general-purpose  digital  computer.  In  addition,  two  Honeywell  124-series 
digital  computers  were  interfaced  with  the  XDS  9300  to  operate  as  parallel 
processors.  This  added  computation  capability  is  needed  for  real-time 
operation  of  the  MIRA  simulation. 

This  section  provides  a brief  introduction  to  the  capabilities  of 
the  MIRA  simulation,  and  documents  its  status  at  the  end  of  this  work 
effort.  In  addition,  the  capability  and  algorithms  added  to  the  simu- 
lation during  the  present  contractual  effort  are  presented  in  detail. 
Included  under  this  topic  are  an  error-navigation  algorithm,  which  iso- 
lates the  effects  of  sensor  errors  on  navigation  performance,  and  auto- 
pilot, altimeter,  and  two-degree-of-f reedom-gyro  models.  A brief  summary 
is  given  of  a transport-aircraft  model  that  has  been  simulated  in  real 
time  to  permit  an  investigation  of  the  MIRA  concept  and  MIRA  configura- 
tions with  this  type  of  vehicle.  Last  of  all  is  a discussion  of  an 
actuator  model  included  in  the  longitudinal  flight-control  system. 
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2.1  MIRA  Simulation  Status 


As  mentioned  previously,  CSDL's  initial  effort  on  the  MIRA  project 
concentrated  upon  the  development  of  a nonreal-time  simulation  capa- 
bility to  allow  the  evaluation  of  the  MIRA  concept  and  MIRA  sensor 
configurations.  A major  goal  of  the  present  effort  was  the  development 
of  an  identical  capability  in  a real-time  digital  simulation.  As  a 
result,  a real-time  simulation  is  available  with  the  following 
features : 

(1)  A tactical  fighter  aircraft  model. 

(2)  Nonlinear  six-degree-of-f reedom  equations  of  motion. 

(3)  A rotating  spherical  earth  model. 

(4)  Nonlinear  aerodynamics  valid  beyond  the  stall  angle-of- 
attack. 

(5)  Three  longitudinal-  and  three  lateral-directional  structural 
modes. 

(6)  A representative  flight-control  system  to  improve  the 
flying  qualities  of  the  vehicle. 

(7)  Turbulence  and  winds. 

(8)  Flight-control  and  navigation-sensor  models.  (These  sensors 
can  be  located  anywhere  on  the  vehicle.  The  available 
navigation-sensor  models  are  laser  gyros  and  quartz-flexure 
accelerometers,  and  can  be  used  for  flight-control  purposes.) 

(9)  An  error-navigation  algorithm. 

(10)  An  instrumented  cockpit. 

Most  of  the  models  included  in  the  simulation  are  discussed  in  detail 

( 1 ) * 

in  the  first-phase  report. 

2.2  Error-Navigation  Algorithm 

The  nonreal-time  version  of  the  MIRA  simulation  contains  a model 

of  a strapdown  local-vertical  wander-azimuth  (LVWA)  navigation  system 

( 2 ) 

obtained  from  the  Standardized  Software  Development  Program. ' With 
this  system,  navigation  performance  is  assessed  by  consideration  of 


Superscript  numerals  refer  to  similarly  numbered  items  in  the  List  of 
References. 
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the  navigation  errors  that  are  obtained  by  differencing  the  outputs 
of  the  navigation-system  model  and  the  corresponding  true  vehicle  states. 
This  scheme  reflects  the  expected  performance  of  an  actual  inertial- 
navigation  system,  since  it  includes  the  effects  of  computational  errors 
in  the  overall  navigation  performance. 

An  alternate  approach  exists  for  evaluating  the  effects  of  the 
flight  environment  on  inertial-sensor  performance,  which  has  been 
incorporated  into  both  the  real-time  and  nonreal-time  versions  of  the 
MIRA  simulation.  It  is  a navigation-  and  attitude-error-propagation 
scheme  using  incremental  velocity-  and  angular-sensor  errors  as  inputs 
together  with  flight-profile  data. ^ In  the  MIPA  simulation,  these 
input  errors  are  obtained  by  differencing  the  output  of  the  simulated 
inertial  sensors  and  the  "ideal"  sensors,  although  inertial-sensor- 
error  models  could  be  used  as  well.  The  outputs  of  the  error-propagation 
algorithm  are  the  errors  in: 

(1)  Position — latitude,  longitude,  altitude. 

(2)  Velocity — north,  east,  vertical. 

(3)  Attitude — roll,  pitch,  and  heading. 

The  main  advantage  resulting  from  the  error-propagation  approach 
is  the  possibility  of  isolating  the  effects  of  the  sensor  errors  on 
navigation  performance.  The  computational  errors  inherent  in  the 
strapdown  LVWA  navigation-system  algorithm  are  eliminated  from  consider- 
ation of  the  navigation  perfromance  when  the  error-navigation  algorithm 
is  employed.  However,  with  both  algorithms  incorporated  into  the 
nonreal-time  MIRA  simulation,  it  is  possible  to  assess  the  magnitude 
of  this  effect. 

The  following  assumptions  were  made  in  the  development  of  the 
navigation-error  equations: 

(1)  The  earth  is  spherical. 

(2)  The  aircraft  equations-of -motion  model  provides  a nominal 
trajectory,  consisting  of  the  body-to-local-vertical  trans- 
formation, the  body  specific  forces,  the  body  angular  rates, 
the  latitude,  longitude,  and  altitude  of  the  vehicle,  and 
the  local  vertical  with  respect  to  earth  velocities. 

(3)  The  inertial-sensor  models  provide  the  errors  in  the  angular 
rates  and  specific  forces. 

Each  of  these  assumptions  is  consistent  with  those  made  in  the  develop- 
ment of  the  MIRA  simulation. 
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2.2.1  Equation  Development 


Since  the  earth  model  used  in  both  the  equations  of  motion  and  the 
strapdown  navigation  model  is  spherical,  a space-stabilized  computational 
frame  is  employed  for  the  navigation-error  equations.  This  has  the  great 
advantage  that  the  navigation  equations,  including  vertical  damping, 
are  perfectly  symmetrical  with  respect  to  each  of  the  three  axes  of  the 
position  and  velocity  variables.  Therefore,  only  one  axis  of  the  three 
need  be  considered  in  developing  the  error  equations,  and  the  results 
extended  to  the  remaining  two  axes  by  symmetry.  The  result  is  that 
the  state  errors  are  arranged  in  a 2 * 3-matrix  form  rather  than  as  a 
6xl-column  vector,  and  the  state  transition  matrix  becomes  2x2  instead 
of  6x6.  Similarly,  the  attitude-error  equations  become  a single-matrix 
vector  equation. 

The  mathematical  flow  diagram  for  the  x-axis  of  the  space-stable 
navigation  algorithm  with  vertical  damping  is  shown  in  Figure  2-1.  The 

corresponding  diagram  for  the  y-axis  is  obtained  by  replacing  AV  , X , 

XI  1 

and  Xj  by  AVy  , Y^ , and  Yj,  respectively,  and  similarly  for  the  z-axis. 


h 


Figure  2-1.  Flow  diagram  for  x-axis  of  space-stable  navigation 
algorithm  with  vertical  damping. 
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The  x-axis  navigation  algorithm  may  be  written  from  Figure  2-1  as 


For  all  three  axes,  the  algorithm  is 

[Rl„  - {[Al  ♦ 1C] chn_x J tR] * [b]  [A»J*  (2-5 


where 


Perturbing  the  elements  of  Eq.  (2-1)  (i.e.,  replacing  X by 

Jn 

(X-  + eX  ),  etc.),  expanding  and  retaining  only  the  first-order  terms 

An  An 
results  in 


Kin  = M Kln-1  +{UA1  + lCl£hn-l(  [Xlln-1  + ™ eAVX_  <2“8> 

*n 


This  is  the  fundamental  error  equation  for  position  and  velocity. 

Each  term  in  Eq.  (2-8). is  now  examined  in  detail.  Consider  first 
the  matrix  [A] , which  may  be  written  as 


• 

1 

> 

ft 

O 

' AtLO 

AtLo' 

[A)  = 

0 

i 

- L 

_AtLO 

AtLO 

2 J 

where 


(2-9) 


(2-10) 


When  altitude  is  perturbed,  L may  be  written  as 


and 


L(h  + eh) 


L(h)  + eL(h) 


= L (h) 


3L 

air 


eh 


(2-11) 


(2-12) 


The  constant  portion  of  the  term  3L/3h  multiplied  by  the  altitude 
error  is  the  altitude-divergence  term  inherent  in  any  undamped  three- 
axis  pure  inertial  navigator.  It  is  the  prime  reason  vertical  damping 
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is  always  incorporated.  The  net  effect  of  altitude  may  be  incorporated 
with  ample  accuracy  by  considering  only  this  first  term,  since  the  rest 
of  the  terms  in  the  series  are  negligible  (e.g.,  1 percent  of  the 
nominal  value  for  altitudes  between  sea  level  and  50,000  feet).  There- 
fore, [A]  may  be  written  as 


[A(h  + eh)] 


[A(h)]  + [EA(h)] 


(2-13) 


where 


[cA(h) ] 


39, 


eh 


At 


LO 


At 


LO 


At 


At 


LO 


LO 


(2-14) 


Next,  the  term  eAVy  has  to  be  expanded  so  it  can  be  evaluated  in 
*1 

terms  of  known  quantities.  It  is  obtained  by  perturbing  the  equation 
for  AVX  . This  latter  quantity  is,  in  turn,  obtained  by  transforming 

the  integrated  outputs  from  the  body  frame  to  the  inertial  frame.  In 
other  words 


[AVt]  = [TB11 [AVj 


(2-15) 


and  perturbing  this  equation  leads  to 


leAVj]  = (Tbi][eAVa]  + [£Tbi][AVa] 


(2-16) 


The  elements  of  the  vector  ( AVft]  are  the  incremental  ideal  vehicle- 
accelerometer  signals,  which  are 


[AVa] 


BI 


fBI 


BI 


(2-17) 


The  ordering  and  signs  of  the  elements  of  [AVA]  are  due  to  the  difference 
in  the  orientation  of  the  axes  between  the  vehicle  and  the  navigation 
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system.  The  vehicle-axes  system  has  a north/east/down  orientation,  while 
the  navigation-axis  system  has  an  up/east/north  orientation. 

The  elements  of  the  vector  [eAVAJ  are  the  accelerometer-error 
signals.  More  specifically,  they  are  the  differences  between  the  incre- 
mental outputs  of  the  actual  accelerometer  models  and  the  ideal  accelerom- 
eter models.  In  particular 


(2-18) 


teGxl  is  a skew-symmetric  matrix  corresponding  to  the  vector  of  inte- 
grated rate  errors,  i.e. 

Eri 

ecJ, 

£P, 

The  elements  of  [cG]  are  obtained  from  the  difference  in  the  angular 
increments  between  the  modeled  gyros  and  the  ideal  gyros  every  cycle 
according  to  the  formulae 

EPcn  ‘ 1 + " PGlJ  + (rdn  * EqGn  > “ (<IGI  * ErG„  > 

n n- 1 n n n n—  1 n n-1 


EqG  * eqG  . + (A0qGn  " qGIn>  - (rGI  ' ePGn  > + ‘PGI  * ErG„  > 
n n-i  n n n n— l n n— 1 


CrG  " erGn  , + (A9rGn  * rGI>  + (qGI  ’ EPGn  > “ <PGI  * EqG„  > 
n n-l  n n n n-l  n n-1 

(2-21) 
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Once  again,  the  differences  in  signs  and  subscripts  are  due  to  the 
different  orientations  of  the  vehicle-  and  navigation-axes  system. 

Equation  (2-16)  may  now  be  rewritten  as 

[eAVj]  = [Tbi1  {[£AVa]  + [eG]  x [AVft]| 


(2-22) 


where  x is  the  cross-product  operation.  Combining  Eq.  (2-22)  and  (2-8), 
and  expanding  to  account  for  the  three  axes,  results  in  the  position- 
and  velocity-vector-error  equation 


[ERln  = IAHeR]^  + {[eA]  + tclehn_i}(R)n-l 


+ tb]  IiTciJITbc1  [IeAVA]  + leG1  x lAVA] 
Use  has  been  made  of  the  equation 


]jJT  (2-23] 


[TBI>  * [TClltTBC] 


(2-24) 


The  major  remaining  step  required  for  the  mechanization  of  Eq.  (2-23) 
is  to  define  [R]  and  [TBJ]  in  terms  of  quantities  available  from  the 
vehicle  model.  In  order  to  do  this,  the  transformation  from  the  local- 
vertical/north  frame  to  the  inertial  frame  is  needed.  This  transformation 
is  given  by 

lTCIJ  * lTAlTfTXJ 


cos 

A 

-sin 

A 

0 * 

'cos 

X 

0 

-sin 

x' 

- 

sin 

A 

cos 

A 

0 

0 

1 

0 

. 0 

0 

1 . 

.sin 

X 

0 

cos 

X. 

Since 

A = wEt  + a 

and  X and  a are  available  from  the  vehicle  model,  each  element  of  this 
transformation  is  defined. 
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Position  components  in  the  inertial  frame  are  obtained  by  trans- 
forming the  radius  vector  from  the  earth  center  to  the  vehicle  center 
of  gravity  (eg)  through  [TCI ] . That  is 


xl‘ 

' R0  + h 

*1 

[TCIJ 

0 

z_ 

0 

I 

m 

(Rq  + h)  cos  X cos  A 
(Rfl  + h)  cos  X cos  A 
(rq  + h)  sin  X 


(2-26) 


Note  that  the  x,y,z-axes  of  the  inertial  frame  are  along  the  earth's 
polar  axis,  positive  north,  along  the  intersection  of  the  initial 
meridian  plane  and  the  equatorial  plane,  positive  out,  and  90  degrees 
east  in  the  equatorial  plane,  respectively.  The  corresponding  axes 
of  the  local-vertical/north  frame  are  up/east/north. 

The  velocity  components  of  the  vehicle  in  the  inertial  frame 
will  now  be  determined.  The  earth-relative  velocities  of  the  vehicle 
in  the  local-vertical/north  frame  are  available  from  the  MIRA  simulation 
and  designated  as  [VCE) . The  velocity  of  the  earth  with  respect  to  the 
inertial  frame  in  the  local-vertical/north  frame  is  designated  lVEI) , 
where 


[VEI]  = [0  u>e(R0  + h)  cos  X 0]T  (2-27) 

The  vehicle  velocity  relative  to  the  inertial  frame  in  local-vertical/ 
north  coordinates  is 


[VCI]  - IVCE]  + ^EI1  (2-28) 

Therefore,  the  vehicle  velocity  in  the  inertial  frame  IVj]  is 

[vll  * (Tci^CI1  (2-29) 

The  elements  of  the  matrix  [R]  are  now  available  from  Eq.  (2-26)  and 
(2-29) . 

The  next  task  is  the  generation  of  the  body-to-inertial  transforma- 
tion, [TB11,  required  to  rotate  the  increment-velocity-error  vector 
into  the  inertial  frame.  ITBI1  may  be  expressed  as  the  product  of  two 
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rotations:  from  the  body  to  the  local-vertical/north  frame,  [ ] ; and 
from  the  local-vertical/north  to  the  inertial  frame,  tTCI).  That  is 


lTBI>  = 


(TClJ  tTBC] 


(2-30) 


A modified  version  of  [T*^]  is  available  from  the  vehicle  equations- 
of-motion  routine.  The  difference  between  what  is  available  and  what 
is  needed  is  due  to  the  difference  in  the  definition  of  the  axes  system 
between  the  body  frame  and  the  navigation  frame.  The  body  frame  is 
defined  in  a forward/right/down  sense,  while  the  navigation  frame  is 
defined  in  an  up/east/north  sense.  Taking  this  difference  in  axes 


^c1 


-T 


the 

elements  of 

[T*c]  are 

BC33 

-T 

BC32 

~T|!C31 

BC23 

T 

BC22 

TbC21 

(2 

BC13 

T 

BC12 

Tbcii 

2.2.2  Attitude  Errors 

Attitude  errors  arise  from  two  different,  but  not  independent, 
sources.  First,  any  errors  in  initial  alignment  or  gyro-drift  compensa- 
tion will  produce  angular  errors  in  the  knowledge  of  the  motion  of  the 
body  frame.  This  information  is  contained  in  [eG] . Second,  accelerom- 
eter and  gyro  compensation  errors,  initial-alignment  errors,  and  mechani- 
zation errors  integrate  into  latitude  and  longitude  errors.  These  lati- 
tude and  longitude  errors  are  angular  errors  in  the  knowledge  of  the 
local-vertical/north  frame,  since  the  computed  attitude  matrix  is  a func- 
tion of  both  where  the  computer  "thinks"  the  body  frame  is,  [eG],  and 
where  it  "thinks"  the  local-vertical/north  frame  is,  [ g A ] . 


The  attitude-error  equation  will  now  be  derived.  First,  the  posi- 
tion errors  must  be  interpreted  as  errors  in  the  transformation  from  the 
inertial  frame  to  the  computed  local-vertical/north  frame.  The  local- 
vertical/north-to-inertial  transformation  was  previously  defined  in 
Eq.  (2-25),  and  its  transpose  is  given  by 


itic'  ■ IVT'ta1 


(2-32) 
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Assume  errors  in  latitude  and  longitudes,  eX  and  eo,  respectively. 

Also,  assume  no  error  in  o)_t  so  that  eo  = eA.  Then  the  computed 

“ a 

transformation  becomes  [Tj^]  where  C is  the  computed  local-vertical/ 
north  frame.  (Tj£]  is  defined  by 

[TIC]  = {[I0]  + eA  [ATx,  (°)  ] t}  [tx1  t 

• { [IQ]  + eo[ATA(0) ]} [Ta]  (2 


(2-33) 


(AT^(O)]1  is  defined  as 


T 3 tTx] 

(AT.  (0)]  = lim  — 

A X--0  aA 


0 0 1 

0 0 0 

-10  0 


(2-34) 


and  (ATa(0) 1 is  defined  as 


(ATa(0)] 


3(Ta1 


0 10 

■10  0 
0 0 0 


(2-35) 


Expanding  Eq.  (2-33),  retaining  only  first-order  terms,  and  noting  that 
[ATx(0)]T  and  [TX)T,  and  [ ATA (0 ) ] and  [TA1  commute  results  in 

(TjA)  * |[IQ]  + eX(ATx(0)]T  + eo(Tx]T(ATA(0) ][TX)| [TIC]  (2-36) 

Postmultiplying  both  sides  of  Eq.  (2-36)  by  (TCI)  and  expanding  leads  to 


[Tc£]  = Ilg]  + “to  cos 


EO  cos  X 


-eo  sin  X 


eo  sin  X (2-37) 


(I0J  + ("Xl 


(2-38) 
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1 
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The  matrix  [ ttx  1 is  the  skew-symmetric  form  of  the  vector 

t it ] = [-eo  sin  X ca  -ea  cos  X]T  (2-39) 

The  total  error  in  the  computed  attitude  matrix  [T££]  may  now  be 
determined.  The  matrix  [Tg£]  is  given  by 

1TBC]  = [Tc£]  [TBcnTBB]  <2-40) 

= {(I0]  + [irX]}lTBC]{(I0]  + [eGX]}  (2-41) 

or 

[crBC]  = " ^BC1 

= [TTX]  [Tbc]  + (Tbc1  [eGx!  (2-42) 

Equation  (2-42)  is  the  matrix  expression  for  the  errors  in  the  computed 
attitude  matrix.  Another  expression  for  the  same  errors  will  be 
derived  as  a function  of  the  attitude  (Euler)  angles,  the  two  equated, 
and  the  attitude  errors  solved  for. 

The  true  attitude  matrix  or  body-to-local-vertical/north  trans- 
formation is  given  by 


[TBC]  = t^l(T9JT[T$)T 


r l 

0 

0 

tv  = 

0 

COS  l() 

sin 

♦ 

.0 

-sin  iji 

cos 

“cos 

6 0 

-sin 

e" 

IT0I  = 

0 

1 

0 

sin 

e o 

cos 

e_ 

“cos 

$ sin 

4> 

o' 

[T*]  = 

-sin 

4 1 cos 

♦ 

0 

0 0 1 


(2-43) 


(2-44) 


(2-45) 


(2-46) 
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If  all  of  the  attitude  angles  are  perturbed  by  small  angles  eip,  cd, 
and  e<j>,  the  resultant  will  be  the  computed  attitude  matrix,  [Tgg] . 
This  matrix  is  given  by 


IT**] 


tVjlIo1  + eiMAT^tO)]}  [T0]T|[Io]  + e6[AT0(O)]T} 

• [T(J,]T|tI0]  + e<t>[ATj)(0)]T}  (2-47) 


[ ATQ (0) ]T  and  [ AT^ (0) ] are  given  by  Eq.  (2-34)  and  (2-35),  respectively, 
while 


[AT^ (0) ] 


lim 


~w~ 


0 0 0 

0 0 1 

.0  -1  o. 


(2-48) 


Expanding  Eq.  (2-47),  retaining  first-order  terms,  and  manipulating 
leads  to 

[eTBC]  = IT£A]  - [Tbc1 

= [TBC]jc<mV0)lT  + (T^]  {eg  [AT0  (0)  )T 

+ [T0  J [AT^O)  HTe]T}[T())]T|  (2-49) 

Equating  £q.  (2-42)  and  (2-49)  gives  the  matrix  form  of  the  attitude- 
error  equation 

e<J>  t AT^  (0)  ] T + [T^]  je6  [AT0  (0)  ]T 

+ eiMTg]  [ AT^  (0 ) ] [T0]T|(T(()]T  = leGxl  + (TC0]  [rrXl  [TBC] 

(2-50) 


Equation  (2-50)  is  equivalent  to  the  vector  equation 

[eG]  + (Tcb](tt]  = e*[e3]  + [T^){e0[e2l  - e«HT0]  [e^}  (2-51) 
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The  right  side  of  Eq.  (2-51)  may  be  expanded  so  that 


[eG]  + [Tcb1  [it  ] 


-cos  0 cos  <li  sin  <(>  0 

cos  9 sin  ij)  cos  cj>  0 

(_  -sin  6 0 1 


' E^  ' 

£0 

.E<t>  . 

(2-52) 


which  can  be  solved  for  the  three  Euler-angle-error  terms  by  inversion 
of  the  matrix.  The  result  is 


■e*‘ 

-cos  <p 

sin  $ 

0 ' 

£ 0 

1 

cos  0 

sin  (J>  cos 

0 

cos  4>  cos 

0 

0 

,E<t>  . 

-cos  ()>  sin 

0 

sin  4>  sin 

0 

-cos  0. 

{(EG]  + tTCB)[7r]} 


(2-53) 


2.2.3  Position  Errors 

Position  errors  in  the  local-vervical/north  frame  for  output  or 
attitude-error  computation  may  be  determined  by  perturbing  the  trigono- 
metric functions  of  a and  X,  and  substituting  the  inertial-position 
components  into  the  trigonometric  functions. 

Latitude  error  is  found  from 


Since 


sin  (X  + eX)  * sin  X + eX  cos  X 


(2-54) 


sin  X 


cos  X 


RE  " R0  + h 
Ah  *=  eRjj 


(2-55) 

(2-56) 

(2-57) 

(2-58) 
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then 


e A * 


eZj  - Ah  sin  A 
R~~ cos- A 


(2-59) 


Similarly,  the  longitude  error  (ec  * eA)  is  found  from 


sin  (A  + eA)  = sin  A + eA  cos  A 


(2-60) 


Since 


sin  A 


K * v? 


(2-61) 


it  follows  that 


eYj  cos  A - eXj  sin  A 
rZ- cos  A 


(2-62) 


XIeYI  - VXI 


(Rg  cos  A)' 


(2-63) 


The  altitude  error  may  be  found  from 


XjEXj  + YjE Yj  + ZjEZj 


(2-64) 


2.2.4  Ground-Velocity  Errors 


The  errors  in  the  local-vertical/north  ground  velocities  are 
obtained  from  Eq.  (2-27),  (2-28),  and  (2-29).  Solving  for  [VCE1  result.B 
in 


IVCE]  - [TIC][V1]  - IVEI] 


(2-65) 


Perturbing  the  parameters  of  this  equation  and  solving  for  the  error 
in  ground  velocities  leads  to 


[eVCE]  = ITIC1  I eVI  1 + UT^HVj]  - [eVEI]  (2-66) 

where  (TIC1  is  obtained  from  Eq.  (2-25);  IeTic)  is  equal  to  Ittx  1 lTIC)  » 
obtained  from  Eq.  (2-37);  [eVj]  is  obtained  from  Eq.  (2-23);  [VjJ  is 
obtained  from  Eq.  (2-29);  and 

IeV__J  = [0  -eXR„  sin  X + eh  cos  X 0]T  (2-67) 

EI  E 

2.2.5  Computational  Frequencies 

The  equations  that  have  been  mechanized  in  the  MIRA  simulation  are 
the  inertial-state-error  equation  (2-23) , the  attitude-error  equation 
(2-53),  the  position-error  equations  (2-59),  (2-63),  (2-64),  and  the 
ground-velocity-error  equation  (2-66) . The  fundamental  computation 
frequency  at  which  these  equations  are  solved  is  5 Hz.  The  only 
operations  that  must  be  performed  at  a higher  frequency  are  those 
required  to  accumulate  the  total  body  misalignment  and  the  incremental 
velocity  errors  in  the  inertial  frame.  This  is  the 

[TBC)j(eAVA)  + [eG]  x [AVa]| 

term  of  Eq.  (2-23) . These  computations  are  performed  at  50  Hz  in  the 
MIRA  simulation.  Ten  samples  of  the  outputs  of  this  computation  are 
summed  and  used  in  the  other  calculations  of  Eq.  (2-23).  In  addition, 
the  midcycle  value  of  (T^j]  is  used  in  the  solution  of  this  equation 
to  obtain  a better  estimate  of  its  average  value  over  a computational 
cycle. 

2 . 3 Autopilot  Algorithm 

Previously,  the  only  possible  inputs  to  the  MIRA  simulation 
vehicle  have  been  the  pilot  commands.  Profile  time  histories  of  the 
desired  signals  are  entered  as  a series  of  points  in  time  along  with 
the  corresponding  amplitudes  of  the  input  variable  at  the  time  points. 

The  computer  program  linearly  interpolates  the  amplitude  of  the  input 
variable  between  the  time  points  to  achieve  the  desired  profiles.  With 
this  capability,  it  is  possible  to  approximate  pilot  inputs  present 
in  flight. 
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However,  with  this  scheme  it  is  difficult  to  have  the  vehicle  fly 
a desired  trajectory  since  it  is  not  always  possible  to  accurately 
determine  the  pilot  commands  that  will  generate  the  trajectory.  To 
overcome  this  shortcoming  in  the  MIRA  simulation,  an  autopilot  has 
been  added  to  allow  the  vehicle  to  fly  a desired  trajectory.  The  auto- 
pilot generates  the  pilot  commands  that  result  in  the  vehicle  following 
the  flight-path,  pitch-attitude-angle,  roll-angle,  heading,  sideslip, 
or  velocity  commands  entered  into  the  simulation.  The  profile-generation 
scheme,  mentioned  in  the  previous  paragraph,  can  be  used  to  generate  the 
desired  trajectory  commands.  The  pilot  commands  that  cause  the  vehicle 
to  fly  the  desired  trajectory  are  generated  based  upon  the  error 
between  the  desired  command  signal  and  the  actual  signal.  In  order 
to  null  out  any  steady-state  errors  that  may  arise,  the  pilot  command 
is  also  a function  of  the  integral  of  the  error  between  the  desired 
and  actual  signal. 

More  specifically,  the  autopilot  equations  programmed  into  the 
MIRA  simulation  are 


s 


x auto 


= K60  69 


+ k/6 e /6e 


+ K6Y 


K/6Y  /6^ 


(2-68) 


where 


66  = ec  - e 

6y  = Yc  - Y 


sy  auto 


K6$  6*  + K/6*  /6* 


(2-69) 


where 

64>  ■ 4>c  — 4> 


6RP  auto  - K6*  + K/6*  /6*  + K6B  66  + K/68  /6B 

where 

6^  ” 4»c  — 4» 

60  - 6C  - 0 


(2-70) 


20 


T auto 


K6V  5V  + K/6V  /6V 


(2-71) 


where 


VC  - VA 


2.4  Altimeter  Model 


Both  the  navigation  algorithm  and  error-navigation  algorithm  incor- 
porated in  the  MIRA  simulation  require  altitude  as  an  input  to  provide 
damping  of  the  vertical  channel.  It  is  possible  to  select  one  of  three 
signals  for  this  use.  Figure  2-2  indicates  the  options.  One  option  is 
true  altitude,  wnich  assumes  a perfect  altimeter.  The  second  option 
results  in  an  altitude  signal,  which  is  generated  by  adding  a bias  to 
the  true  altitude  and  lagging  the  sum.  The  last  option  requires 
atmospheric  pressure,  which  is  generated  in  the  MIRA  simulation  by 


pA(ib/ft2) 


CONVERSION 

lb/ft2  -» in.  Hg 
0.014139 


ALTITUDE  GENERATION 
EQUATIONS 

(Eq.  (2-72),  (2-73),  or  (2-741) 


Figure  2-2.  Altimeter  model. 
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table  lookup  as  a function  of  true  altitude.  The  atmospheric  signal 

is  lagged  and  a bias  added  to  simulate  the  operation  of  the  sensor. 

The  air-data  computer  is  simulated  by  generating  altitude  from  atmos- 

(4) 

pheric  pressure  using  the  following  functions.  ' 

For  h < 36089.23  (6.68323876  < p.) 

p — — A 

1,  p v 0.190263105 ) 

1 \29. 9212598/  j (2-72) 

For  36089.23  < h_  < 65616.79  (1.61673069  < P.  < 6.68323876) 

— p — A 

hp  = 36089.23  + 20805.8257  in  |iiM|23876j  (2-73) 

For  65616.79  < h < 104986.87  (p.  < 1.61673069) 

— p — A 

|p  .-0.029271247) 

1 -(l. 61673069/  j <2"74> 

2.5  Stabilator-Actuator  Model 

The  flight-control-system  model  incorporated  into  the  MIRA  simula- 
tion was  developed  for  the  Survivable  Flight  Control  System  (SFCS) 
Program. This  system  was  selected  because  it  is  consistent  with 
the  structural-mode  data  used  in  the  simulation,  and  because  it  has 
been  verified  through  flight  test.  The  system  incorporated  into  the 
MIRA  simulation  was  simplified  somewhat  in  that  a complex  stabilator- 
actuator  model  was  replaced  by  a simpler  second-order- system  model. 

During  the  present  simulation  development  effort,  the  SFCS  actuator 
model,  shown  in  Figure  2-3,  has  been  incorporated  into  the  simulation. 
The  actuator  model  shown  differs  from  that  shown  in  Reference  5,  due 
to  the  elimination  of  some  very  high  frequency  filters,  which  have  no 
effect  on  the  MIRA  simulation. 

2.6  Two-Degree-of-Freedom-Gyro  Model 

Models  of  laser  gyros  were  incorporated  into  the  MIRA  simulation 
during  the  initial  phase  of  the  program.  This  capability  has  been 
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expanded  to  include  models  of  two-degree-of-f reedom  gyroscopes.  Three 
models,  representing  a channel  of  a two-degree-of-f reedom  gyro  instru- 
ment mounted  along  each  axis  of  a strapdown  orthogonal  triad  parallel 
to  the  vehicle-body  axes,  are  now  included  in  the  MIRA  simulation.  A 
block  diagram  of  one  channel  of  a two-degree-of-f reedom  gyro  is  shown 
in  Figure  2-4.  The  process  begins  with  the  limiting  and  integration 
of  the  roll-rate  signal  of  the  vehicle,  including  bending  effects. 

The  coupling  errors  due  to  the  misalignment  of  the  gyro  are  then  added 
in  and  the  resulting  quantity  converted  to  pulses  by  multiplying  by 
the  scale  factor. 

Additional  error  terms  are  summed  with  the  resulting  signal  at 

2 

this  junction.  The  error  effects  included  are:  a bias  error,  g,  g , 
gxg-  and  uxu-dependent  errors.  All  possible  sources  of  error  are 
included  in  the  two-degree-of-f reedom-gyro  model  to  allow  the  greatest 
flexibility  in  exercising  the  MIRA  simulation.  For  example,  a term 
may  have  an  effect  for  a particular  orientation  of  a two-degree-of- 
freedom  gyro,  while  having  a negligible  effect  for  a different  orienta- 
tion. 

The  two-degree-of-freedom-gyro  signal  (including  the  error  effects) 
is  then  lagged  and  quantized  with  the  integer  part  converted  from  pulses  to 
radians  via  the  scale  factor.  The  scale-factor  error  is  incorporated  at 
this  time.  The  output  of  the  gyro  model  is  then  available  for  use  in  other 
subroutines  of  the  MIRA  simulation,  e.g.,  the  navigation-system  model, 
the  flight-control  system,  and  the  error-navigation  algorithm. 

2.7  Transport  Model 

A model  of  a transport-type  aircraft  has  been  simulated  in  real 
time  to  allow  its  investigation  with  the  MIRA  concept  and  MIRA  configura- 
tions. The  aerodynamics  of  this  vehicle  and  flight-control  system 
included  in  the  simulation  are  documented  in  detail  in  Reference  7 . The 
details  of  these  subjects  will  not  be  discussed  in  this  report  for 
the  sake  of  brevity  and  economy.  However,  a general  description  of  the 
transport  model  is  presented  to  give  the  reader  a brief  idea  of  the 
features  that  it  incorporates. 

The  aerodynamics  of  the  transport  vehicle  modeled  are  valid  over 
the  whole  flight  regime  of  the  aircraft  from  landing  approach  to  cruise. 


Figure  2-4.  Roll  channel  of  two-degree-of-freedom-gyro  model 


A table-lookup  and  linear-interpolation  scheme  is  used  to  obtain  the 
proper  aerodynamic  forces  and  moments  of  the  vehicle.  The  transport 
model  includes  a number  of  control  surfaces: 

(1)  A horizontal  stabilizer  for  pitch  trim. 

(2)  An  elevator  for  dynamic  pitch  control. 

(3)  Inboard  and  outboard  ailerons  and  spoilers  for  roll  control. 

(4)  A rudder  for  yaw  control. 

(5)  Speed  brakes,  flaps,  and  landing  gtar. 

An  augmentation  system  has  also  been  incorporated  into  the  simulation. 
It  includes  a yaw  damper  to  increase  the  damping  ratio  of  the  Dutch- 
roll  mode.  The  transport  simulation  has  been  configured  to  operate  in 
conjunction  with  the  transport-vehicle  cockpit  available  at  CSDL  to 
give  added  realism  to  the  environment  in  which  the  pilot  flies  the 
simulated  vehicle. 
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SECTION  3 

MIRA  CONFIGURATION  EVALUATION 


This  section  is  devoted  to  an  evaluation  of  potential  MIRA  con- 
figurations using  the  simulation  capability  outlined  in  Section  2.1. 

The  data  runs  were  selected  to  provide  a maximum  amount  of  information 
with  a minimum  number  of  runs.  From  the  data  collected,  it  is  possible 
to  assess  the  comparative  effects  of: 

(1)  A realistic  and  ideal  environment. 

(2)  Lever-arm  effects  due  to  sensor  location. 

(3)  Two-degree-of-f reedom  and  laser  gyros. 

(4)  Different  sensor  configurations. 

(5)  The  effect  of  the  sensor  configurations  on  navigation  and 
flight  control. 

(6)  The  use  of  navigation  sensors  for  flight-control  purposes. 

The  initial  part  of  this  section  is  devoted  to  a discussion  of  the 
simulation  input  data,  i.e.,  sensor  configurations,  flight  trajectories, 
sensor-error  coefficients,  structural-mode  effects  on  the  sensors,  etc. 
The  results  obtained  via  the  MIRA  simulation  are  then  presented, 
followed  by  a discussion  of  the  simulation  results  and  conclusions 
drawn  from  the  results. 

3. 1 Simulation  Input  Data 

Two  distinct  sets  of  simulation  runs  were  performed  under  this 
effort  as  outlined  in  Table  3-1.  One  set  of  runs  is  short-term 
(30-second)  and  was  designed  to  excite  the  longitudinal  and  lateral- 
directional  transient  modes  of  the  vehicle,  so  that  the  comparative 
effects  of  the  potential  MIRA  configurations  on  the  handling  qualities 
of  the  vehicle  could  be  investigated.  A second  set  of  runs  of  1-hour 
. duration  was  designed  to  evaluate  the  navigation  performance  of  the 

potential  MIRA  configurations. 
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Table  3-1.  MIRA  configuration  evaluation  runs. 


Configuration 

Gyro  Type 

Sensor 

Location 

Turbulence 

and 

Modes 

Short 

Run 

Long 

Run 

Navigation 

Sensors 

for 

Flight 

Control 

Triad 

Laser 

FS313 

No 

Yes 

Yes 

No 

Triad 

Two-degree-of- 

freedom 

FS313 

No 

Yes 

Yes 

No 

Triad 

Laser 

FS313 

Yes 

Yes 

Yes 

No 

Triad 

Laser 

FS77 

Yes 

Yes 

Yes 

No 

Triad 

Two-degree-of- 

freedom 

FS313 

Yes 

Yes 

Yes 

No 

Triad 

Two-degree-of- 

freedom 

FS77 

Yes 

Yes 

Yes 

No 

Dodecahedron 

Laser 

FS313 

Yes 

Yes 

Yes 

No 

Octahedron 

Two-degree-of- 

freedom 

FS313 

Yes 

Yes 

Yes 

No 

Triad 

Laser 

FS313 

Yes 

Yes 

No 

Yes 

Triad 

Two-degree-of- 

freedom 

FS313 

Yes 

Yes 

No 

Yes 

The  nominal  flight  condition  of  the  vehicle  for  all  of  the  simu- 
lation runs  is  0.5  Mach  at  an  altitude  of  5000  feet.  Figure  3-1  depicts 
the  autopilot-trajectory  commands  for  the  short-duration  simulation 
runs.  The  bank-angle  and  heading  commands  were  selected  to  cause  the 
vehicle  to  undergo  a coordinated  30-degree  banked  turn  until  the  heading 
angle  changed  by  30  degrees.  The  5-degree  change  in  flight-path  angle 
was  arbitrarily  selected  to  excite  the  longitudinal  modes  of  the  vehicle. 
For  the  long-duration  navigation  evaluation  runs,  the  autopilot  commands 
were  selected  so  that  the  vehicle  would  fly  due  north  at  a constant 
speed  and  constant  altitude.  When  turbulence  is  present,  its  nominal 
rms  value  is  6 feet  per  second,  which  corresponds  to  the  level  of  clear- 
air  turbulence  specified  in  Reference  8 at  the  nominal  altitude. 

The  orthogonal  triad  was  selected  as  one  of  the  basic  sensor  con- 
figurations for  evaluation  in  this  study.  Results  were  obtained  for 
this  sensor  configuration  using  either  laser  or  two-degree-of-freedom 
gyros,  located  near  the  eg  of  the  vehicle  (FS313)  or  just  in  front  of 
the  pilot's  station  (FS77) . In  addition,  the  comparative  effects  of 
the  environment  were  evaluated  using  this  configuration  by  performing 
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Figure  3-1. 


Trajectory  commands  for  short-duration 
evaluation  runs. 


the  simulation  runs  in  a turbulent  environment  with  the  structural 
modes  excited,  and  in  a nonturbulent  environment  solving  only  the  rigid- 
body  equations  of  the  vehicle.  The  effect  of  using  navigation  sensors 
for  flight-control  purposes  has  also  been  assessed  using  the  orthogonal- 
triad  configuration  with  laser  or  two-degree-of-f reedom  gyros. 

Another  sensor  configuration  selected  for  evaluation  was  the 
Strapdown  Inertial  Reference  Unit  (SIRU) , which  employs  six  single- 
degree-of-f reedom  gyros  and  six  linear  accelerometers.  The  instrument 
input  axes  are  arranged  in  a nonorthogonal  pattern  that  corresponds  to 
the  array  of  normals  to  the  faces  of  a regular  dodecahedron  as  shown 
in  Figure  3-2.  This  array  is  completely  symmetrical,  with  the  acute 
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Figure  3-2.  The  dodecahedron  configuration. 

angle  between  any  two  axes  equal  to  63.4  degrees.  Since  all  input  axes 
are  nonorthogonal , each  instrument's  output  contains  a measure  of  redun- 
dant data. 

The  question  arises  as  to  how  a skewed  sensor  configuration  is 
evaluated  using  the  MIRA  simulation,  which  has  only  a laser  or  two- 
degree-of-freedom  gyro  and  accelerometer  triad  modeled.  The  approach 
selected  uses  the  covariance  of  the  errors  of  the  sensor  outputs,  which 

are  computed  relative  to  an  orthogonal  triad.  For  example,  the  error- 

2 

covariance  matrix  of  an  orthogonal  triad  is  a [Iq] , assuming  identical 

instruments  along  each  axis.  It  is  shown  in  Reference  9 that,  relative 

to  an  orthogonal  triad,  the  error-covariance  matrix  of  the  dodecahedron 

2 

configuration  is  0.5oe[IQ].  A nominal  set  of  error  coefficients  has 
been  selected  for  the  instruments  of  the  orthogonal  triad.  These 
coefficients  are  adjusted  by  the  square  root  of  0.5  to  account  for  the 
effect  due  to  the  orientation  of  the  sensors  in  the  dodecahedron  con- 
figuration. 

The  last  configuration  selected  for  evaluation  uses  four  two- 
degree-of-freedom  gyros  with  the  gyro  spin  axes  each  normal  to  one  of 
the  faces  of  a regular  semi-octahedron.  The  sensor-measurement  axes 
are  oriented  such  that  they  all  lie  on  the  surface  of  a cone  where  the 
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the  error-covariance  matrix  is  0.612ae[I0J,  which  is  quite  close  to 
that  of  the  dodecahedron  configuration. Therefore,  when  evaluating 
this  configuration  with  the  MIRA  simulation,  the  two-degree-of- freedom 
error  coefficients  of  the  orthogonal  triad  were  adjusted  by  the  square 
root  of  0.612. 

In  the  evaluation  of  the  dodecahedron-  and  octahedron-sensor  con- 
figurations, the  accelerometer-  and  gyro-error  coefficients  have  been 
adjusted  by  the  same  gain  factor.  This  implies  that  the  gyros  and 
accelerometers  are  mounted  in  the  same  fashion  (e.g.,  eight  linear 
accelerometers  are  employed  in  the  octahedron  configuration) . However, 
the  results  obtained  for  this  particular  case  are  applicable  to  other 
configurations  as  well.  Consider,  for  example,  the  configuration  shown 
in  Figure  3-5,  which  uses  only  five  accelerometers.  The  input  axes  of 
the  instruments  are  distributed  on  a cone,  each  axis  separated  from  its 
neighbor  by  72  degrees  in  azimuth  about  the  cone  axis.  The  cone  axis 
is  the  roll  axis,  and  lies  at  an  angle  of  54  degrees,  44  minutes,  8 sec- 
onds from  the  surface  of  the  cone.  The  error-covariance  matrix  for 
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Figure  3-5.  Conical  configuration. 


32 


2 (7i 

this  configuration  is  0.6oe[Ig],  which  is  extremely  close  to  that 
obtained  for  the  octahedron  configuration.  Thus,  the  result  obtained 
from  the  MIRA  simulation  for  the  octahedron  configuration  would  be 
applicable  to  the  conical  configuration  as  well. 

All  of  the  configurations  selected  were  evaluated  with  the  sensors 
colocated  at  FS313.  This  location  is  just  forward  of  the  vehicle  eg, 
which  is  at  FS320  for  the  particular  vehicle  configuration  used. 

Table  3-2  gives  the  distances  along  the  x-  and  z-body  axes  from  the  eg 
to  the  sensor  package  for  this  sensor  location.  This  table  also  con- 
tains the  body-axes  distances  from  the  vehicle  eg  to  FS77  where  the 
orthogonal  triad  configuration  containing  laser  or  two-degree-of- 
freedom  gyros  was  evaluated  to  obtain  an  estimate  of  lever-arm  effects 
upon  the  potential  MIRA  configurations.  The  pilot's  station  is  located 
at  FS119  as  a point  of  reference.  Both  of  these  sensor  locations  were 
recommended  in  Reference  8 based  upon  the  availability  of  space  and 
the  structural-mode  characteristics  of  the  vehicle. 

Table  3-2.  Sensor  location  relative  to 
vehicle  center  of  gravity. 


The  effects  of  the  three  longitudinal-  and  three  lateral-directional 
structural  modes  on  the  sensor  responses  are  given  in  Table  3-3  for  the 
two  sensor  locations  just  discussed.  Parameters  are  given  for  the  nor- 
mal and  lateral  accelerometers,  and  the  roll-,  pitch-,  and  yaw-rate 
gyros.  This  information  was  taken  from  References  5 and  6,  which  do  not 
include  the  effects  of  the  modes  on  the  longitudinal  accelerometer  since 
this  instrument  is  not  used  for  flight-control  purposes.  Therefore, 
the  effect  of  the  inodes  on  this  sensor  has  been  assumed  to  be  negligible. 
The  effect  of  the  modes  on  the  gyro  signals  is  given  by  the  product  of 
the  coefficients  in  Table  3-3  and  the  mode  rate  in  radians  per  second, 
while  the  effect  on  the  accelerometer  signals  is  given  by  the  product 
of  the  appropriate  coefficients  and  the  mode  acceleration  in  rad/s^. 

The  nominal  frequencies  of  the  longitudinal  modes  are  71.0,  87.5,  and 
141.0  radians  per  second,  while  those  of  the  lateral-directional  modes 
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Table  3-3.  Aeroelastic  coefficients  of  sensor  responses 


Longitudinal  Sensors 

Coefficient 

Coefficient 

Coefficient 

Sensor 

Location 

Units 

of  First 
Longitudinal 

of  Second 
Longitudinal 

of  Third 
Longitudinal 

Mode 

Mode 

Mode 

Normal 

FS77 

ft 

4. 774  (10~3) 

-1. 467 (lo"3) 

-8. 333 (IQ-4) 

Accelerometer 

FS313 

ft 

-1.840(lo'3) 

5. 006 (lo"3) 

-3.075  (lO-3) 

Pitch-Rate 

FS77 

rad 

-3.0608(10'4) 

2. 239 (lo"3) 

-1.14 (lo'4) 

Gyro 

FS313 

rad 

-2.33  <10'4) 

2.64  <10"4) 

0.0 

Lateral- 

Directional  Sensors 

Coefficient 

Coefficient 

Coefficient 

of  First 

of  Second 

of  Third 

Sensor 

Location 

Units 

Lateral- 

Lateral- 

Lateral- 

Directional 

Directional 

Directional 

Mode 

Mode 

Mode 

Lateral 

FS77 

ft 

0.0 

0.0 

-3.885  do'4) 

Accelerometer 

FS313 

ft 

0.0 

0.0 

9. 066  do'5) 

Roll-Rate 

FS77 

rad 

-5.5  (lo'4) 

-4 

6.3  (10  ) 

2.9  dO'3) 

Gyro 

FS313 

rad 

-3.  33  (lo"4) 

5.  333  (lO"4) 

2.267(1 0'3) 

Yaw- Rate 

FS77 

0.0 

0.0 

-1 . 4 ( 10~3) 

Gyro 

FS313 

HI 

0.0 

0.0 

-3. 5 (10~4) 

are  39.7,  67.6,  and  80.8  radians  per  second.  The  lateral  accelerometer 
and  yaw-rate  gyro  are  affected  by  only  the  highest  frequency  lateral- 
directional  mode,  while  the  roll-rate-gyro  signal  is  influenced  by  all 
three  lateral-directional  modes.  The  three  longitudinal-structural 
modes  modeled  influence  both  the  pitch-rate  gyro  and  normal  accelerom- 
eter signals. 

A nominal  set  of  parameters  was  selected  for  the  laser-gyro, 
quartz-flexure-accelerometer,  and  two-degree-of-f reedom-gyro  models 
for  the  MIRA  evaluation  runs.  These  parameters  are  given  in  Tables  3-4 
3-5,  and  3-6.  The  nominal  parameter  set  was  used  for  each  instrument 
of  the  orthogonal  triad  programmed  in  the  MIRA  simulation,  with  the 
sign  of  the  error  parameters  selected  to  be  positive  or  negative  with 
equal  probability  to  introduce  some  randomness  into  the  parameter- 
selection  process.  The  same  set  of  accelerometer  parameters  was  used 
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Table  3-4.  Ring- laser-gvro  nominal  parameters 


f 

I 


I 

> 


[ 


L 


Parameter 

Value 

Scale  Factor 

Bias 

Scale-Factor  Error 

Misalignment  Coefficients 

1.57  sec/pulse 

0.01  deg/h 

5 ppm 

5 (10-5)  rad 

Table  3-5.  Quartz-flexure-accelerometer 
nominal  parameters. 


Parameter 

Value 

Scale  Factor 

32,000  pps/g 

Bias 

50  yg 

Scale-Factor  Error 

100  ppm 

Misalignment  Coefficients 

5 (10'S)  rad 

Cross-Coupling 

.5 

Coefficient 

1.96(10  ) rad/g 

Scale-Factor  Nonlinearity 

Coefficient 

40  yg/g 

with  either  the  laser  or  two-degree-of-f reedom  gyros.  For  configurations 
other  than  the  orthogonal  triad,  the  magnitude  of  the  error  coefficients 
was  adjusted  by  the  square  root  of  the  coefficient  of  the  error-covariance 
matrix. 

It  has  been  assumed  for  the  orthogonal  triad  with  two-degree-of- 
freedom  gryos  that  one  instrument  is  aligned  with  the  x-instrument  out- 
put axis  along  the  x-body  axis,  to  measure  the  roll  rate;  and  the  y-body 
axis,  to  measure  the  pitch-rate  signal.  The  spin  axis  of  this  instru- 
ment is  down,  parallel  to  the  vehicle  z-body  axis.  The  yaw-rate  signal 
was  obtained  from  the  y-output  axis  of  a two-degree-of-freedom  gyro 
mounted  with  the  x-instrument  axis  parallel  to  the  y-body  axis,  the 
y-instrument  axis  parallel  to  the  z-body  axis,  and  the  gyro  spin  axis 
parallel  to  the  x-body  axis.  The  x-axis  output  signal  of  this  gyro  was 
not  used  in  the  simulation.  It  has  been  assumed  that  the  instruments 
are  mounted  on  the  vehicle  in  this  fashion  to  minimize  the  effects  of 
the  g-dependent  errors  on  the  gyro  output  signals. 
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Table  3-6.  Two-degree-of-f reedom-gvro  nominal  parameters 


Parameter 

Units 

Scale  Factor 

sec/pulse 

Bias 

deg/h 

Scale-Factor 

Error 

ppm 

Misalignment 

Coefficients 

rad 

g-De pendent 
Errors 

deg/h/g 

2 

g -Dependent 
Errors 

deg/h/g2 

gxg-Dependent 

Errors 

deg/h/g2 

u>xw-  Dependent 
Errors 

deg/h/ (rad/s) 2 

x-Output-Axis 

Value 


y-Output-Axis 

Value 


Source 


“V 

for  x-output 

"V  “gz 

for  y-output 


9y 

9z 


"x  9z 


“g*'  “gy 
“gy'  “g2 


3.2  Simulation  Results 

This  section  contains  the  results  obtained  for  the  cases  des- 
cribed in  Section  3.1,  using  the  CSDL-designed  and  developed  MIRA 
simulation.  The  1-hour  runs  designed  to  evaluate  the  navigation  per- 
formance of  the  candidate  MIRA  systems  will  be  discussed  first. 

Figure  3-6  shows  the  system  responses  obtained  with  the  orthogonal 
triad  configuration  using  two-degree-of-f reedom  gyros  located  near  the 
eg  of  the  vehicle  at  FS313.  The  results  obtained  for  this  case  are 
typical  of  those  obtained  for  the  other  cases  in  the  study. 

The  autopilot  causes  the  vehicle  to  follow  the  desired  true- 
north  trajectory  at  constant  speed  and  altitude  very  well  in  the 
turbulent  environment  over  the  duration  of  the  flight.  The  Euler 
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MIRA  simulation  responses:  triad  configuration  at  FS313; 
two-degree-of-f reedom  gyros;  1-hour  run  in  turbulent 
environment  (sheet  1 of  8) . 
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Figure  3-6.  MIRA  simulation  responses:  triad  configuration  at  F3313 
two-degree-of- freedom  gyros;  1-hour  run  in  turbulent 
environment  (sheet  2 of  8) . 
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Figure  3-6.  rtIRA  simulation  responses.  triad  configuration  at  FS313 
two-degree-of- freedom  gyros;  1-hour  run  in  turbulent 
environment  (sheet  5 of  8) . 
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Figure  3-6.  MIRA  simulation  responses:  triad  configuration  at  FS313; 

two-degree-of-freedom  gyros;  1-hour  run  in  turbulent 
environment  (sheet  6 of  8) . 
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Figure  3-6.  MIRA  simulation  responses:  triad  configuration  at  FS313, 
two-degree-of- freedom  gyros;  1-hour  run  in  turbulent 
environment  (sheet  7 of  8) . 
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Figure  3-6.  MIRA  simulation  responses:  triad  configuration  at  F3313 
two-degree-of-f reedom  gyros:  1-hour  run  in  turbulent 
environment  (sheet  8 of  8) . 
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angles  deviate  less  than  ±1.5  degrees  from  their  trim  values,  as  do 
the  angle  of  attack  and  sideslip  angle.  The  speed  of  the  vehicle  is 
maintained  within  1 foot  per  second,  while  the  maximum  deviation  of 
the  altitude  from  the  5000-foot  trim  value  is  17  feet.  Other  flight- 
profile  variables  of  interest  are  the  flight-path  angle,  which  is  less 
than  1 degree  in  magnitude,  and  the  altitude  rate  and  east  ground  vel- 
ocity, which  are  less  than  9 feet  per  second  in  magnitude.  The  hori- 
zontal tail  and  and  throttle  deflections,  commanded  to  keep  the  vehicle 
on  its  trajectory,  are  of  relatively  small  magnitude.  The  commanded 
lateral-directional  controls  are  of  much  larger  magnitude,  approximately 
5 degrees  maximum  for  the  rudder,  and  2 degrees  maximum  for  the  aileron. 

However,  these  are  not  unreasonable  values.  Also  presented  in  Figure  3-6 
are  the  flight-control-sensor  time  histories.  The  yaw-  and  roll-rate 
signals  are  less  than  2.5  degrees  per  second,  and  the  pitch-rate-gyro 
signal  less  than  1 degree  per  second.  The  incremental  normal  acceler- 
ation is  on  the  order  of  ±0.2  g,  while  the  lateral  acceleration  peaks 
at  0.05  g. 

The  structural-mode  accelerations,  rates,  and  shapes  are  included 
in  Figure  3-6  to  give  an  indication  of  the  magnitude  and  character  of 
the  excitations.  The  level  of  mode  excitation  is  very  nearly  the  same 
for  all  of  the  simulation  runs.  The  effect  of  the  modes  on  the  sensor 
signals  can  be  estimated  from  these  time  histories  and  the  data  con- 
tained in  Table  3-3.  For  example,  the  effect  of  both  the  first  and 
second  longitudinal  modes  on  the  normal  accelerometer  signal  is  less 
than  0.006  g.  The  effect  of  the  third  longitudinal  mode  is  even  less, 
about  0.001  g.  The  maximum  contribution  to  the  pitch-rate-gyro  signal 
is  caused  by  the  first  longitudinal  mode,  and  its  magnitude  is  0.024 
degree  per  second.  The  roll-rate  gyro  is  affected  almost  equally  by 
the  three  lateral-directional  modes  with  the  respective  peak  contribu- 
tion to  the  signal  0.1,  0.06,  0.1  degree  per  second.  The  yaw-rate  gyro 
is  affected  by  the  third  lateral-directional  mode  only,  and  the  peak 
level  of  its  effect  is  less  than  0.01  degree  per  second.  The  same  is 
true  of  the  lateral  accelerometer  with  the  magnitude  of  the  effect  less 
than  0.0005  g. 

Figure  3-6  also  includes  the  output  responses  of  the  navigation- 
system  model  included  in  the  MIRA  simulation.  These  responses  appear 
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to  be  almost  duplicates  of  their  vehicle  counterparts.  Perhaps  the 
most  noticeable  differences  are  evident  in  the  east  ground  velocity, 
which  has  a long-term  bias  in  the  navigation-system  output;  and  longi- 
tude error,  which  is  four  times  larger  in  magnitude,  although  it  is 
an  extremely  small  quantity. 

The  capability  and  accuracy  of  the  navigation  system  is  better 
understood  and  evaluated  by  considering  error  responses  along  with  the 
outputs  of  the  system  rather  than  just  the  system  responses  alone.  With 
this  in  mind,  two  sets  of  navigation-error  variables  have  been  included 
in  Figure  3-6.  One  set  is  obtained  by  differencing  the  navigation- 
system  outputs  and  the  vehicle  responses.  The  corresponding  set  of 
error  signals  obtained  from  the  error-navigation  algorithm  (and  due 
solely  to  the  sensor  errors)  is  also  included  for  comparison.  On  the 
whole,  the  two  sets  of  error  responses  agree  remarkably  well  for  this 
case,  being  almost  identical  overlays  in  most  cases.  The  altitude- 
and  altitude-rate-error  signals  do  not  agree,  but  this  is  due  to  the 
fact  that  the  navigation  system  includes  a third-order  loop  in  the 
vertical  channel,  while  that  in  the  error-navigation  algorithm  is  only 
second  order.  The  navigation-system  errors  tend  to  be  noisier  due  to 
the  differencing  with  the  true  state  of  the  vehicle.  Although  not  evi- 
dent from  the  responses  of  Figure  3-6,  the  Schuler  period  for  the  sys- 
tem shown  is  of  the  correct  magnitude  as  determined  from  a 2-hour  rather 
than  a 1-hour  run. 

The  comparative  navigation  performance  of  the  other  MIRA  systems 
investigated  was  determined  from  the  system-error  responses  obtained  by 
differencing  the  navigation-system  and  vehicle  outputs,  as  well  as 
from  the  outputs  of  the  error-navigation  algorithm.  These  responses 
are  presented  in  Figures  3-7  through  3-13  for  the  configuration  pre- 
sented previously  in  Table  3-1.  The  simulation  runs  in  these  figures 
are  again  1-hour  long. 

The  error  responses  given  in  Figure  3-7  represent  a baseline 
case  for  the  simulation  runs  performed  in  that  they  were  obtained  with 
an  orthogonal-triad  configuration  located  at  FS313  in  a nonturbulent 
environment  with  only  the  vehicle  rigid-body  equations  being  solved. 

The  structural-mode  equations  of  the  vehicle  were  not  solved.  The 
responses  for  the  corresponding  system  in  a turbulent  environment 
with  the  structural  modes  of  the  vehicle  being  excited  are  given  in 
Figure  3-8. 
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Figure  3-7.  MIRA  simulation  responses:  triad  configuration  at  FS313 
laser  gyros;  1-hour  run  in  nonturbulent  environment 
without  structural  modes  (sheet  2 of  2) . 
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Figure  3-8.  MIRA  simulation  responses:  triad  configuration  at  FS313; 
laser  gyros;  1-hour  run  in  turbulent  environment 
(sheet  1 of  2)  . 
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MIRA  simulation  responses:  triad  configuration  at  FS313 
two-degree-of- freedom  gyros;  1-hour  run  in  nonturbulent 
environment  without  structural  modes  (sheet  1 of  2) . 
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Figure  3-«.  MIRA  simulation  responses:  triad  configuration  at  FS313; 

two-degree-of-freedom  gyros;  1-hour  run  in  nonturbulent 
environment  without  structural  modes  (sheet  2 of  2) . 
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MIRA  simulation  responses:  triad  configuration  at  FS77 
laser  gyros;  1-hour  run  in  turbulent  environment 
(sheet  1 of  2) . 
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MIRA  simulation  responses:  triad  configuration  at  FS77 
laser  gyros;  1-hour  run  in  turbulent  environment 
(sheet  2 of  2) . 
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MIRA  simulation  responses:  triad  configuration  at  FS77 
two-degree-of-freedom  gyros;  1-hour  run  in  turbulent 
environment  (sheet  1 of  2) . 
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Figure  3-11.  MIRA  simulation  responses:  triad  configuration  at  FS77; 

two-degree-of-freedom  gyros;  1-hour  run  in  turbulent 
environment  (sheet  2 of  2) . 
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MIRA  simulation  responses:  dodecahedron  configuration 
at  FS313;  laser  gyros;  1-hour  run  in  turbulent 
environment  (sheet  1 of  2) . 
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MIRA  simulation  responses:  dodecahedron  confi* 
at  FS313;  laser  gyros;  1-hour  run  in  turbulent 
environment  (sheet  2 of  2) . 
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MIRA  simulation  responses:  octahedron  configuration  at 
FS313;  two-degree-of-freedom  gyros;  1-hour  run  in 
turbulent  environment  (sheet  1 of  2)  . 
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Figure  3-13 


MIRA  simulation  responses:  octahedron  configuration  at 
FS313;  two-degree-of- freedom  gyros;  1-hour  run  in 
turbulent  environment  (sheet  2 of  2) . 


A comparison  of  the  responses  in  these  two  figures  leads  to  the 
conclusion  that  the  error-navigation  outputs  obtained  are  almost  exact 
duplicates  of  each  other.  That  is,  the  navigation  errors  due  to  the 
sensor  errors  are  the  same  in  both  environments.  The  maximum  error  in 
latitude  and  longitude  due  to  the  sensor  errors  is  about  0.4  nautical 
mile.  The  errors  obtained  by  differencing  the  navigation-system  and 
vehicle  responses  are  the  same  as  the  responses  of  the  error-navigation 
algorithm  outputs  for  the  benign-environment  case.  This  verifies  the 
operation  and  development  of  the  navigation-system  algorithm.  The 
responses  obtained  for  the  turbulent  environment  are  slightly  noisier, 
as  would  be  expected.  The  effect  of  the  turbulence  and  structural 
modes  on  the  navigation  accuracy  can  be  determined  by  comparing  the 
navigation- system  errors  and  error-navigation  responses  of  Figure  3-8. 
First  of  all,  the  navigation-system  Euler  angle  errors  tend  to  be  quite 
a bit  noisier  due  to  the  high  frequencies  present  in  the  vehicle  Euler 
angles.  Perhaps  the  most  dramatic  difference  between  the  two  sets  of 
error  variables  is  in  the  latitude  response.  This  variable  is  biased 
negatively  for  the  navigation-system  error  rather  than  oscillating  about 
zero,  as  in  the  case  of  the  error-navigation-algorithm  output.  A maxi- 
mum latitude  error  of  0.7  nautical  mile  results.  A bias  is  present  in 
the  north  ground-velocity  navigation-system  error  as  well. 

Using  Figures  3-9  and  3-6,  a similar  comparison  may  be  made  for 
an  orthogonal-triad-sensor  configuration  comprised  of  two-degree-of- 
freedom  gyros.  The  responses  of  Figure  3-9  are  for  the  benign  environ- 
ment, while  the  responses  of  Figure  3-6  have  already  been  discussed. 

In  the  nonturbulent  environment,  the  navigation-system  errors  and  the 
error-algorithm  outputs  match  extremely  well.  The  sensors  introduce 
latitude  errors  on  the  order  of  0.9  nautical  mile,  while  the  longitude 
error  is  one-half  of  that.  The  errors  obtained  with  the  two-degree- 
of-freedom-gyro  system  appear  to  be  slightly  larger  than  those  obtained 
with  the  laser-gyro  triad  presented  in  Figure  3-7.  More  specifically, 
most  of  the  differences  appear  to  be  due  to  different  drift  rates 
associated  with  the  sensor  models  that  result  in  the  laser  and  two- 
degree-of-freedom  triad  responses  being  skewed  with  respect  to  each 
other. 

A comparison  of  the  responses  obtained  with  the  two-degree-of- 
freedom-gyro  triads  in  the  calm  and  turbulent  environment  leads  to  the 
conclusion  that  the  error-navigation-algorithm  outputs  are  basically  the 
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same  in  each  case.  The  error-navigation  outputs  in  the  turbulent  en- 
vironment tend  to  be  noisier,  particularly  the  Euler-angle  errors,  but 
the  mean  signal  tend  to  be  the  same.  The  navigation  system  errors  in 
the  turbulent  environment  are  noisier  yet  due  to  the  high  frequencies 
present  in  the  vehicle  variables.  Some  of  the  navigation-system  errors 
are  also  larger,  such  as  latitude  and  north  ground  velocity. 

A comparison  of  the  responses  in  the  turbulent  environment  for 
the  laser-gyro  and  two-degree-of-f reedom-gyro  triad  configurations 
(Figures  3-8  and  3-6)  reveals  that  the  two-degree-of-freedom-gyro 
model  introduces  much  more  noise  into  the  output  signal,  particularly 
the  Euler-angle-error  signals.  In  addition,  the  latitude  error  is  on 
the  order  of  three  times  larger  with  the  two-degree-of-freedom-gyro 
sensors.  It  is  on  the  order  of  2.5  nautical  miles  for  the  configura- 
tion chosen. 

The  lever-arm  effects  on  navigation  performance  may  be  assessed 
by  comparing  Figures  3-10  and  3-11  with  the  figures  just  discussed. 

Figure  3-10  shows  the  error  responses  obtained  at  FS77  with  the  laser- 
gyro  orthogonal  triad,  while  Figure  3-11  shows  the  corresponding 
responses  for  the  orthogonal  triad  with  two-degree-of-f reedom  gyros. 

As  for  the  sensor-error-induced  navigation  errors,  they  tend  to  be  the 
same  without  respect  to  sensor  location  for  the  trajectory  flown.  This 
is  true  for  both  the  two-degree-of-freedom-  and  laser-gyro  configura- 
tions. As  far  as  the  navigation  system  is  concerned,  it  tends  to  intro- 
duce noise  into  the  system  as  evidenced  by  the  noisier  error-navigation- 
algorithm  responses.  This  is  particularly  true  of  the  ground-velocity- 
error  signals.  In  addition,  the  latitude  error  is  larger  as  a result 
of  the  navigation  algorithm,  and  a bias  is  introduced  into  the  north- 
ground-velocity  error. 

The  last  two  sets  of  navigation-error  responses  obtained  were  those 
for  the  dodecahedron  and  octahedron  sensor  configurations  located  at 
FS313.  Figures  3-12  and  3-13  show  these  results.  A comparison  of  the 
dodecahedron-configuration  responses  with  the  orthogonal  laser-triad 
responses  reveals  that  the  error  responses  for  the  dodecahedron  are 
essentially  scaled-down  versions  of  their  laser-triad  counterparts. 

This  is  as  it  should  be  since  the  error  coefficients  of  the  sensor 
models  for  the  dodecahedron  configurations  were  scaled  down  by  the 
square  root  of  the  error  covariance  gain.  The  only  noticeably  differ- 
ent response  is  the  error-navigation-algorithm  attitude  error,  which 
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has  a different  sign  for  the  two  cases  of  interest.  The  same  conclu- 
sions are  reached  when  the  results  for  the  octahedron  configuration 
are  compared  to  the  laser-triad  results.  A direct  comparison  of  the 
results  obtained  for  the  dodecahedron  and  octahedron  configurations 
reveals  that  the  dodecahedron  responses  are  cleaner  and  smaller  than 
those  of  the  octahedron. 

A further  comparison  of  the  navigation  performance  of  the  various 
sensor  configurations  considered  is  provided  in  Tables  3-7,  3-8,  and 
3-9.  These  tables  contain  the  maximum  magnitude,  mean,  and  standard 
deviation  of  the  error  signals  obtained  for  the  1-hour  simulation 
runs.  A consideration  of  the  maximum  magnitudes  of  the  various  error 
signals  leads  to  the  following  conclusions: 

(1)  The  values  of  the  error-navigation-algorithm  maximums  are 
independent  of  environment  and  position  for  the  laser-triad 
configuration. 

(2)  The  maximum  values  of  the  error-navigation  outputs  for  the 
two-degree-of-freedom-gyro  triad  tend  to  be  the  same  for 
both  FS313  and  FS77.  However,  these  quantities  tend  to  be 
larger  than  those  obtained  in  the  benign  environment. 

(3)  The  maximum  values  of  the  navigation-system-  and  error- 
algorithm  outputs  are  the  same  in  the  calm  environment. 

(4)  The  smallest  maximum  errors  are  obtained  for  the  dodecahedron 
configuration.  The  laser  orthogonal  triad  appears  to  follow, 
being  slightly  better  than  the  octahedron  configuration. 

(5)  The  octahedron  configuration  produces  smaller  maximum  errors 
than  the  triad  with  two-degree-of ^freedom  gyros. 

(6)  On  the  average,  in  a turbulent  environment  the  maximum  mag- 
nitudes of  the  navigation-system  errors  tend  to  be  about 
50-percent  larger  than  their  error-navigation-algorithm 
counterparts. 

The  results  available  from  Table  3-8  (the  listing  of  the  mean 
navigation  errors)  tend  to  be  less  conclusive  and  fewer  in  number  than 
those  obtained  by  considering  the  magnitude  of  the  maximum  errors. 
However,  some  of  the  same  conclusions  are  evident.  For  example, 

(1)  The  mean  errors  obtained  with  the  laser-triad  configuration 
appear  to  be  independent  of  position  and  environment. 
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Table  3-7.  The  maximum  magnitude  of  the  navigation-error  signals  for  the  1-hour  simulation  runs. 


Table  3-8.  The  mean  of  navigation-error  signals  for  the  1-hour  simulation  runs. 


■0.4427  -0.4345  0.24097  0.1255  -0.3130  0.18843  -0.5384  0.1223 


Table  3-9.  The  standard  deviation  of  the  navigation-error  signals  for  the  1-hour  simulation  runs. 


0.7944  0.7913  0.8136  0.81727  0.4326  0.4953  0.7817  0.7869 


(2)  The  smallest  mean  errors  result  with  the  dodecahedron  con- 
figuration. However,  the  next  best  configuration  is  not 
clear. 

(3)  The  navigation-system  errors  and  the  error-navigation- 
algorithm  outputs  have  the  same  mean  value  in  the  nonturbu- 
lent  environment. 

(4)  The  mean  errors  obtained  with  the  octahedron  configuration 
are  smaller  than  those  obtained  with  the  two-degree-of- 
freedom-triad  configuration. 

(5)  In  the  calm  environment,  the  mean  values  of  the  laser-triad- 
error  variables  are  smaller  than  those  of  the  triad  comprised 
of  two-degree-of- freedom  gyros. 

(6)  The  mean  errors  obtained  with  the  two-degree-of- freedom-gyro 
triad  are  a function  of  position  and  environment. 

(7)  The  largest  mean  errors  are  obtained  with  the  sensors  located 
at  FS313  in  a turbulent  environment,  while  the  smallest 
errors  are  obtained  at  FS313  without  the  turbulence. 

Table  3-9  (the  standard  deviations  of  the  error  variables) 
provides  many  conclusions  concerning  the  relative  merits  of  the 
navigation  performance  of  the  sensor  configurations  evaluated.  The 
following  conclusions  may  be  drawn  from  the  table: 

(1)  The  standard  deviations  of  the  error  variables  associated 
with  the  laser  configurations  are  less  than  the  corresponding 
entries  for  the  two-degree-of-freedom  gyros. 

(2)  The  quantities  obtained  with  the  error-navigation  algorithm 
for  the  laser-gyro  configuration  do  not  appear  to  be 
affected  greatly  by  position  and  environment.  However, 
when  the  standard  deviations  of  the  navigation  output  errors 
are  considered,  this  does  not  appear  to  be  true  to  the 

same  extent. 

(3)  The  error-algorithm  outputs  obtained  for  the  two-degree-of- 
freedom-triad  configuration  show  a relatively  large  varia- 
tion with  position  and  environment. 

(4)  The  dodecahedron  configuration  appears  to  produce  the  smallest 
standard  deviations.  The  next  best  configuration,  when  all 

of  the  results  are  taken  into  consideration,  is  not  clear. 
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(5)  The  octahedron  configuration  produced  better  results  than  the 
two-degree-of-freedom  orthogonal  triad. 

The  results  of  the  short-duration  simulation  runs  will  now  be 
considered.  As  mentioned  previously,  the  commanded  trajectory  was 
selected  to  excite  the  transient  modes  of  the  system.  This  was  done 
by  having  the  vehicle  undergo  a 30-degree  banked  turn  until  a 30-degree 
change  in  heading  was  achieved.  The  wings  were  then  leveled  and  a 
5-degree  change  in  flight-path  angle  commanded. 

Figure  3-14  contains  a complete  set  of  time  histories  obtained 
for  one  of  the  short-duration  runs.  The  flight-control  and  navigation 
sensors  are  all  colocated  at  FS313.  The  autopilot  causes  the  vehicle 
to  follow  the  commanded  trajectory  reasonably  well.  The  bank  angle 
settles  at  30  degrees  after  a slight  initial  overshoot,  and  the  heading 
increases  linearly  until  a 30-degree  change  is  achieved.  The  speed  is 
maintained  within  3 feet  per  second  of  the  trim  value,  and  it  deviates 
by  about  that  much  only  when  the  change  in  flight-path  angle  is  commanded. 
The  flight-path  angle  achieves  a peak  of  8 degrees,  and  this  is  due  to 
the  coupling  in  the  vehicle  responses.  The  leveling  of  the  vehicle's 
wings  after  the  desired  heading  is  achieved  causes  a transient  in  the 
flight-path  angle.  The  5-degree  commanded  change  in  flight-path  angle 
comes  when  this  transient  is  occurring.  The  turn  is  fairly  well  coordi- 
nated with  a side  acceleration  of  less  than  0.1  g.  A normal  accelera- 
tion of  0.8  g results  when  the  pull-up  command  is  initiated.  The 
flight-control-sensor  signals  are  oscillatory,  reflecting  the  vehicle 
modal  responses  in  the  turbulent  environment.  The  vehicle-control 
deflections  are  reasonable  with  a 12.0-degree  aileron  pulse  required 
to  achieve  the  30-degree  bank  angle,  while  the  rudder  deflection  is 
essentially  parabolic  in  shape,  reflecting  a command  to  maintain  a 
constant  yaw  rate  and  zero  sideslip  angle.  The  longitudinal  controls 
are  excited  primarily  when  the  flight-path-angle  command  is  initiated. 

More  thrust  is  called  for  to  maintain  speed,  while  the  horizontal  tail 
deflection  causes  the  vehicle  to  pitch  up. 

The  time  histories  obtained  for  the  calm  environment,  which 
correspond  to  those  just  discussed,  have  also  been  obtained.  They  are 
just  smoothed  versions  of  the  responses  of  Figure  3-14  and  are  not 
presented  for  this  reason. 

The  structural-mode  responses  are  also  included  in  Figure  3-14 
to  provide  an  indication  of  the  magnitude  and  character  of  their 
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Figure  3-14 


MIRA  simulation  responses:  all  sensors  located  at  FS313 
laser  navigation  gyros;  30-second  run  in  turbulent 
environment  (sheet  1 of  4) . 
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MIRA  simulation  responses:  all  sensors  located  at  FS313 
laser  navigation  gyros;  30-second  run  in  turbulent 
environment  (sheet  2 of  4)  . 
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excitation.  The  mode  excitation  results  in  the  following  maximum 
effect  or  the  sensors: 

(1)  0.004  g on  the  normal  accelerometer  and  0.002  degree  per 
second  on  the  pitch-rate  gyro,  each  due  to  the  first 
longitudinal  mode. 

(2)  4.0(10  5)  g on  the  lateral  accelerometer  and  0.007  degree 
per  second  on  the  yaw-rate-gyro  signal,  each  due  to  the 
third  lateral-directional  mode. 

(3)  0.1  degree  per  second  on  the  roll-rate  gyro,  due  to  the 
first  lateral-directional  mode. 

None  of  these  effects  is  of  major  importance. 

The  last  set  of  variables  included  in  Figure  3-14  are  the  outputs 
of  the  navigation-system  model  and  the  error  responses  due  to  the 
differencing  of  the  navigation  system,  and  vehicle  responses  and  the 
output  of  the  error-navigation  algorithm.  First,  the  error-navigation 
outputs  are  almost  identical  to  those  obtained  in  the  nonturbulent 
environment  without  structural  modes.  It  is  also  evident  that  the 
environment  and  dynamic  behavior  of  the  vehicle  cause  large  navigation- 
system  errors  relative  to  those  caused  by  the  sensor  errors. 

The  MIRA  simulation  responses  included  in  Figure  3-15  correspond 
to  those  of  Figure  3-14  with  the  exception  that  the  sensors  are  located 
at  FS77  rather  than  near  the  eg  of  the  vehicle.  There  are  several 
differences  evident  from  these  two  sets  of  responses.  One  of  them  is 
that  the  damping  of  the  lateral-directional  responses  is  greater  when 
the  sensors  are  located  at  FS77,  resulting  in  less  oscillatory 
repponses.  On  the  other  hand,  the  longitudinal  responses  are  about 
the  same  in  this  regard.  The  lateral-directional  structural  modes 
are  also  excited  to  a slightly  greater  degree  when  the  sensors  are 
located  near  the  pilot's  station.  However,  this  change  is  not  signifi- 
cant enought  to  cause  a major  effect  on  the  sensors.  A comparison  of 
Figures  3-14  and  3-15  shows  that  the  navigation  errors  due  to  the 
sensor  errors  are  independent  of  sensor  location.  The  navigation-system 
errors  are  larger  when  the  sensors  are  located  at  FS77,  reflecting  the 
lever-arm  effects  present  in  this  case. 

The  next  two  sets  of  time  histories  presented  are  the  same  as 
those  just  discussed  with  the  exception  that  two-degree-of-freedom- 
gyro  models  were  used  instead  of  laser-gyro  models.  The  responses  in 
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MIRA  simulation  responses:  all  sensors  located  at  FS77 
laser  navigation  gyros;  30-second  run  in  turbulent 
environment  (sheet  1 of  4) . 


Figure  3-15 


MIRA  simulation  responses:  all  sensors  located  at  FS77 
laser  navigation  gyros;  30-second  run  in  turbulent 
environment  (sheet  2 of  4) . 
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Figure  3-15.  MIRA  simulation  responses:  all  sensors  located  at  FS77 
laser  navigation  gyros;  30-second  run  in  turbulent 
environment  (sheet  3 of  4). 
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Figure  3-16  were  obtained  with  the  sensors  located  at  FS313,  while  those 
of  Figure  3-17  were  obtained  with  the  sensors  at  FS77.  Since  the  flight- 
control-sensor  models  were  used  in  all  cases,  the  vehicle  responses, 
structural-mode  responses,  etc.  are  identical  for  the  corresponding  cases. 
Hence,  they  are  not  repeated  in  Figures  3-16  and  3-17;  only  the  naviga- 
tion errors  are  presented.  Comparing  the  responses  obtained  with  the 
sensors  located  near  the  eg  reveals  that  the  sensor-generated  navigator 
errors  are  larger  with  the  two-degree-of-freedom-gyro  models.  The  Euler- 
angle  and  longitude  errors  are  an  order  of  magnitude  larger.  When  the 
navigation-system  errors  are  examined,  the  differences  in  the  Euler- 
angle  errors  stand  out.  With  the  two-degree~of-freedom  gyros,  the 
navigation-system  errors  and  error-navigation-algorithm  responses  are 
identical  for  the  Euler  angles,  and  are  much  larger  than  those  obtained 
with  the  laser-gyro  models.  When  Figure  3-17  is  considered  (i.e.,  the 
triad  is  comprised  of  two-degree-of-freedom-gyro  models  located  at 
FS77) , the  error-algorithm  responses  duplicate  those  obtained  at 
FS313,  implying  that  the  sensor-generated  navigation  errors  are  not 
a strong  function  of  position  and  environment.  The  navigation-system 
errors  are,  however,  much  larger,  duplicating  those  of  the  corresponding 
laser-gyro  case  except  for  the  Euler  angles.  This  implies  that  the 
lever-arm  effects  play  a major  part  in  the  generation  of  the  navigation 
errors  compared  to  those  of  the  sensor  errors. 

Figure  3-18  and  3-19  contain  responses  representing  the  incorpora- 
tion of  the  dodecahedron  and  octagon  navigation-sensor  configurations 
into  the  vehicle.  The  basic  vehicle  responses  are  the  same  as  those 
presented  in  Figure  3-14  since  the  flight-control-sensor  models  are 
used  in  both  cases.  The  errors  shown  in  Figures  3-18  and  3-19  are 
just  scaled  versions  of  the  responses  obtained  with  the  triad  configura- 
tion. The  responses  obtained  for  the  dodecahedron  configuration 
correspond  to  those  of  Figure  3-14  (i.e.,  the  triad  with  laser  gyros), 
while  the  responses  obtained  for  the  octahedron  configuration  correspond 
to  those  of  Figure  3-16  obtained  for  an  orthogonal  triad  with  the 
two-degree-of-freedom-gyro  models.  The  relative  improvement  is  greater 
with  the  dodecahedron  configuration,  as  would  be  expected  since  its 
error-covariance  gain  is  less  than  that  obtained  for  the  octahedron 
conf iguration . 

The  last  system  to  be  evaluated  used  navigation  sensors  for 
flight-control  as  well  as  navigation  purposes.  The  resulting  responses 
are  presented  in  Figure  3-20.  These  responses,  obtained  for  an 
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sensor  configuration;  30-second  run 
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Figure  3-20.  MIRA  simulation  responses:  all  sensors  located  at  FS313 
laser  navigation  gyros — navigation  sensors  used  for 
flight  control;  30-second  run  in  turbulent  environment 
(sheet  1 of  4) . 
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orthogonal-triad  configuration  of  a laser-gyro  model  located  at  FS313 , 
correspond  to  those  presented  in  Figure  3-14.  The  differences  among 
the  two  sets  of  responses  are  not  very  dramatic,  and  the  basic  vehicle 
responses  do  not  differ  significantly,  if  at  all.  Some  differences 
are  evident  in  the  flight-control  signals,  witn  those  of  Figure  3-20 
showing  the  effect  of  the  sensor  quantization.  The  structural-mode 
responses  are  essentially  the  same  in  both  cases,  with  those  obtained 
for  the  dual-purpose  navigation  sensors  being  perhaps  a slight  bit 
larger  in  magnitude.  The  maximum  value  of  the  accelerations  is  slightly 
higher,  resulting  from  the  use  of  the  navigation  sensors  for  flight 
control,  but  the  effect  is  certainly  nothing  that  is  going  to  be 
discernible  to  the  pilot.  Similar  conclusions  can  be  drawn  for  the 
orthogonal-triad  configuration  with  the  two-degree-of- freedom-gyro 
model.  However,  the  responses  for  this  case  are  not  shown,  since 
they  add  nothing  to  the  information  already  presented.  Therefore, 
there  does  not  appear  to  be  any  strong  technical  reason  why  the  naviga- 
tion sensors  cannot  be  used  for  flight-control  purposes. 


SECTION  4 


TEST  AND  DEVELOPMENT  PLAN 


This  section  presents  and  discusses  the  topics  that  should  be 
included  in  a test  and  development  plan  for  MIRA- type  sensor  systems. 
Several  of  the  items  included  in  the  plan  are  normally  a part  of  test 
plans  for  conventional  sensor  systems,  thus  forming  a base  upon  which 
to  build  the  test  plan.  However,  the  uniqueness  of  the  MIRA  concept 
dictates  that  a much  more  comprehensive  and  detailed  test  plan  be 
developed  and  carried  out.  It  is  this  aspect  of  the  proposed  test 
plan  that  will  receive  the  major  consideration.  A combination  of 
analyses,  laboratory  and  van  testing,  and  use  of  the  MIRA  simulation 
is  proposed  to  characterize  the  parameters  of  a MIRA  system  and  ensure 
its  proper  operation. 

The  initial  requirement  for  the  test  plan  is  a basic  characteri- 
zation of  the  instruments  included  in  the  MIRA  system.  This  aspect 
is  no  different  from  that  which  would  be  done  in  the  testing  of  any 
conventional  sensor  system,  and  includes  laboratory  tests  to  determine 
the  basic  parameters  of  the  sensors  (e.g.,  scale  factors,  rate  limits, 
etc.).  The  next  item  that  should  be  covered  in  the  test  plan  is  a 
determination  of  the  nominal  navigation  and  flight-control  performance 
achievable  with  the  sensor  system.  This  would  involve  a series  of 
laboratory  and  van  tests  to  ascertain  how  well  the  sensor  system 
meets  the  system  specification  and  requirements  of  the  dynamic  environ- 
ment in  which  the  vehicle  will  be  flying.  The  MIRA  simulation  can  be 
useful  in  this  area  because  the  basic  sensor  characteristics  established 
for  this  testing  can  be  incorporated  into  it,  and  the  navigation  and 
flight-control  performance  achieved  in  the  test  can  be  verified  with 
the  simulation. 

The  next  requirement  that  must  be  determined  during  the  testing 
of  the  MIRA  system  is  a definition  of  the  relative  differences  existing 
among  the  system  sensors.  This  information  is  needed  to  establish  the 
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basic  failure  detection  and  isolation  (FDI)  capability  of  the  MIRA 
system.  The  parameters  that  would  have  to  be  determined  during  this 
phase  of  the  testing  include  misalignment  effects  and  those  effects 
associated  with  a dynamic  environment  (e.g.,  the  magnitude  of  the 
g-dependent  effects).  The  dynamic  response  of  the  individual  sensors 
will  also  be  different  due  to  different  bandwidths,  nonlinearities, 
and  the  like.  The  noise  level  associated  with  the  instrument,  which 
may  be  small  in  a dynamic  environment,  must  also  be  determined  to 
characterize  the  quiescent  level  of  the  instrument ' s operation.  An 
extensive  combination  of  laboratory  and  van  testing  in  conjunction 
with  analyses  will  be -needed  to  determine  the  required  parameters  and 
the  degree  of  confidence  required  for  reliable  operation  of  the  MIRA 
system. 

The  next  logical  step  in  the  test  plan  will  be  a determination  of 
the  failure  modes  of  the  sensor  system.  A combination  of  analysis 
and  appropriate  testing  to  confirm  the  analysis  will  be  required. 

Once  the  failures  have  been  catalogued,  an  investigation  of  the  effect 
of  the  failure  modes  on  the  performance  of  the  MIRA  and  vehicle  systems 
should  be  undertaken.  This  phase  of  the  testing  should  concentrate 
upon  two  aspects  of  the  problem.  One  deals  with  hard  failures  and 
their  effect  upon  the  flight -control  system  and  transient  behavior 
of  the  aircraft  to  ensure  that  the  safety  of  the  vehicle  is  not 
compromised;  the  effect  of  the  failures  on  the  handling  qualities  of 
the  vehicle  is  an  area  of  prime  consideration.  The  other  aspect  of  the 
problem  concerns  the  long-term  navigation  performance  of  the  system 
with  instrument  failures.  The  intent  of  the  tests  performed  at  this 
time  should  be  to  gather  statistical  data  for  identifying  the  effects 
of  soft  failures  leading  to  a set  of  FDI  requirements.  Once  again, 
a combination  of  van  testing  and  verification  by  simulation  should  be 
required  during  this  phase  of  the  program. 

The  next  tests  required  are  those  designed  to  specify  the 
parameters  and  verify  the  operation  of  the  FDI  mechanism  of  the  sensor 
system.  With  the  knowledge  obtained  in  the  previous  tests,  it  should 
be  possible  to  induce  failures  in  the  sensor  system  and  verify  that 
the  failure  is  detected  and  isolated,  and  the  proper  reconfiguration 
is  undertaken.  Laboatory  and  van  testing  will  again  be  required.  A 
tool  such  as  the  MIRA  simulation  could  also  be  used  to  confirm  the 
results  of  the  testing  and  provide  more  exhaustive  and  comprehensive 
verification  of  the  FDI  algorithm. 
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The  last  item  that  should  be  considered  in  the  test  plan  is  the 
effect  of  the  sensor  reconfiguration  on  the  performance  of  the  vehicle. 
A.  particular  example  would  be  the  transient  induced  into  the  control 
system  due  to  the  reconfiguration  of  the  sensors.  A combination  of 
analyses  and  use  of  the  MIRA  simulation  should  provide  answers  to 
many  of  the  questions  that  arise  in  this  regard. 
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SECTION  5 

SUMMARY  AND  CONCLUSIONS 


The  MIRAT  technical  effort  at  CSDL  has  been  dedicated  to  the 
accomplishment  of  three  major  goals: 

(1)  The  continuation  of  the  simulation  development  effort  begun 
under  the  first  CSDL  MIRA  program. 

(2)  The  utilization  of  the  simulation  tools  in  the  test  and 
evaluation  phase  of  the  MIRA  program. 

(3)  The  consideration  of  a test  and  development  plan  for 
MIRA-type  systems. 

< 

The  desired  goals  have  been  achieved  in  each  of  these  three  areas. 

The  simulation  development  phase  of  CSDL's  effort  concentrated 
on  the  conversion  of  the  nonreal-time  simulation  (developed  during 
the  first  phase  of  the  CSDL  MIRA  effort)  to  digital  real-time  operation. 
This  effort  was  given  first  priority  in  the  program  plan,  and  the 
result  is  a simulation  capability  that  can  be  used  to  evaluate  a wide 
class  of  MIRA  systems  in  many  different  environments.  The  simulation 
is  coded  in  fixed-point  assembly  language,  and  has  been  designed  to 
operate  on  an  XDS  9300  general-purpose  digital  computer,  and  two 
Honeywell  124  series  computers,  which  operate  as  parallel  processors. 

The  capability  of  the  MIRA  simulation  has  also  been  increased 
significantly  during  the  course  of  this  program  by  the  addition  of 
several  subroutines.  A navigation  and  attitude  error-propagation 
algorithm  using  the  incremental-velocity  and  angular-sensor  errors  as 
inputs  together  with  flight -profile  data  has  been  incorporated  into  the 
MIRA  simulation.  It  is  now  possible  to  isolate  the  effects  of  sensor 
errors  on  navigation  performance  with  this  capability.  Models  of 
two-degree-of-freedom  gyros  were  also  incorporated  into  the  MIRA 
• simulation  to  augment  the  laser-gyro  models  previously  included.  A 
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model  of  a transport  vehicle  was  developed  to  permit  the  investigation 
of  the  operation  of  sensor  systems  in  vehicles  of  this  type.  Last  of 
all,  an  air-data  sensor  model  was  added  to  support  the  proposed  MIRA 
conf igurations . 

The  simulation  runs  performed  during  this  program  were  designed 
with  two  goals  in  mind.  Parametric  studies  were  performed  to  provide 
an  indication  of  the  effects  of  sensor  type,  sensor  location,  sensor 
errors,  etc.  In  addition,  specific  potential  MIRA  configurations 
were  evaluated  to  provide  an  assessment  of  their  capability.  This 
phase  of  the  effort  was  hampered  by  the  fact  that  CSDL  did  not  receive 
the  MIRA  configurations  selected  by  McDonnell  Douglas  for  evaluation  as 
called  for  in  the  contract  Statement  of  Work.  In  lieu  of  this,  CSDL 
selected  and  evaluated  sensor  configurations  that  were  suitable  for 
MIRA  applications. 

Two  major  sets  of  computer  runs  were  performed.  Eight  1-hour 
runs  were  made  to  determine  the  navigation  performance  of  the  various 
configurations  selected,  and  ten  30-second  runs  were  made  to  determine 
the  effect  of  the  configurations  on  the  transient  behavior  of  the 
aircraft  systems. 

The  specific  configurations  evaluated  were  an  orthogonal  triad 
with  laser  gyros  located  near  the  eg  of  the  vehicle  and  then  forward  of 
the  pilot's  station;  the  same  configuration  with  two-degree-of-freedom- 
gyro  models;  and  dodecahedron  and  octahedron  configurations.  Other 
variables  in  the  simulation  were  the  environment  and  the  use  of  the 
navigation  sensors  for  flight-control  purposes. 

Many  conclusions  were  drawn  from  the  simulation  runs  (see 
Section  3).  The  results  concern  the  effects  of  the  parametric  investi- 
gation and  the  relative  merits  of  the  configurations  evaluated.  How- 
ever, it  is  felt  that  the  simulation  runs  performed  and  potential 
MIRA  configurations  evaluated  were  not  extensive  enough  to  justify  the 
selection  and  recommendation  of  the  most  suitable  specific-sensor 
configuration  for  the  MIRA  application. 

Many  suggestions  were  made  for  items  that  should  be  included  in 
a test  and  development  plan  for  MIRA  systems.  The  basic  plan  calls  for 
a combination  of  analyses,  laboratory  and  van  testing,  and  simulation 
verification  to  ensure  the  proper  operation  of  a MIRA-type  system. 

This  plan  includes  the  characterization  of  sensor  instruments  and  the 
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determination  of  system  performance  and  failure  modes,  operations 
associated  with  the  testing  of  sensor  systems  in  general.  However, 
because  MIRA  systems  require  a much  more  extensive  plan,  suggestions 
were  made  to  include  differential  instrument  measurements,  the  effect 
of  the  failure  modes  on  flight-control  and  navigation  performance,  the 
verification  of  the  FDI  and  reconfiguration  performance  of  the  system, 
and  a determination  of  the  effects  of  the  transients  introduced  into 
the  vehicle  system  due  to  the  sensor  reconfiguration. 
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